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Abstract  This  paper  presents  a  novel  method  for  convert¬ 
ing  any  unstructured  quadrilateral  or  hexahedral  mesh  to  a 
generalized  T-spline  surface  or  solid  T-spline,  based  on  the 
rational  T-spline  basis  functions.  Our  conversion  algorithm 
consists  of  two  stages:  the  topology  stage  and  the  geome¬ 
try  stage.  In  the  topology  stage,  the  input  quadrilateral  or 
hexahedral  mesh  is  taken  as  the  initial  T-mesh.  To  construct 
a  gap-free  T-spline,  templates  are  designed  for  each  type  of 
node  and  applied  to  elements  in  the  input  mesh.  In  the  geom¬ 
etry  stage,  an  efficient  surface  fitting  technique  is  developed 
to  improve  the  surface  accuracy  with  sharp  feature  preser¬ 
vation.  The  constructed  T-spline  surface  and  solid  T-spline 
interpolate  every  boundary  node  in  the  input  mesh,  with  C2- 
continuity  everywhere  except  the  local  region  around  irreg¬ 
ular  nodes.  Finally,  a  Bezier  extraction  technique  is  devel¬ 
oped  and  linear  independence  of  the  constructed  T- splines  is 
studied  to  facilitate  T-spline  based  isogeometric  analysis. 
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1  Introduction 

For  the  sake  of  integration  of  engineering  design  and  analy¬ 
sis,  isogeometric  analysis  was  proposed  [4,  1]  which  utilizes 
NURBS  (Non-Uniform  Rational  B -Spline)  or  T- splines  as  a 
basis.  Over  the  forty-year  history  of  commercial  finite  ele¬ 
ment  analysis,  there  are  a  large  number  of  polygonal  meshes 
accumulated  with  the  development  of  automatic  mesh  gen¬ 
eration  techniques  and  finite  element  analysis  technology 
based  on  polygonal  meshes.  For  example,  Figure  1(a)  shows 
an  unstructured  hexahedral  mesh  of  a  gear  assembly.  Hence, 
a  solution  to  converting  these  polygonal  meshes  to  T-splines 
is  needed,  which  provides  engineers  with  the  opportunity 
to  transition  legacy  bilinear  quadrilateral  surface  meshes  or 
trilinear  hexahedral  meshes  to  T-splines,  analyze  them  us¬ 
ing  isogeometric  analysis  and  compare  the  results  with  tra¬ 
ditional  finite  element  technology.  In  addition,  the  T-spline 
representation  is  a  more  compact  way  to  represent  geometry 
compared  with  polygonal  meshes  and  has  better  continuity. 

Previous  approaches  on  converting  meshes  to  spline  rep¬ 
resentations  involved  approximating  the  data  by  determining 
the  topology  and  choosing  a  parameterization.  In  [5,  6],  a 
conversion  method  from  a  triangle  mesh  of  arbitrary  topol¬ 
ogy  into  a  T-spline  surface  was  proposed  based  on  peri¬ 
odic  global  parameterization.  A  polycube  map,  which  mim¬ 
ics  the  input  mesh  in  a  topologically  correct  and  geometri¬ 
cally  meaningful  manner,  was  utilized  as  parametric  domain 
to  construct  T-splines  [14]. 

In  our  earlier  work,  we  developed  an  algorithm  for  con¬ 
verting  any  unstructured  quadrilateral  mesh  to  a  standard  T- 
spline,  whose  basis  functions  form  a  partition  of  unity  [15]. 
In  this  method,  many  nodes  need  to  be  inserted  in  order  to 
make  the  T-spline  standard.  To  reduce  the  number  of  in¬ 
serted  nodes,  as  a  follow-up  we  generalize  the  T-spline  def¬ 
inition  to  the  rational  T-spline  in  this  paper.  The  new  ra¬ 
tional  T-spline  basis  functions  have  the  property  of  a  par- 
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Fig.  1  The  gear  assembly,  (a)  The  input  unstructured  hexahedral  mesh;  (b)  the  constructed  solid  T-spline  and  T-mesh;  and  (c)  the  extracted  solid 
Bezier  elements  with  some  elements  removed  to  show  the  interior  mesh,  (d-f)  show  details. 


tition  of  unity  not  only  for  standard  T- splines,  but  also  for 
semi-standard  and  non-standard  T-splines.  Here,  we  focus 
on  converting  an  arbitrary  unstructured  quadrilateral  or  hex¬ 
ahedral  mesh  to  a  rational  bicubic  T-spline  surface  or  tricu- 
bic  solid  T-spline.  There  are  two  main  stages  in  the  conver¬ 
sion  algorithm:  the  topology  stage  and  the  geometry  stage. 
We  take  the  input  mesh  directly  as  the  initial  T-mesh,  and 
the  topology  stage  aims  to  make  the  initial  T-mesh  gap-free 
by  designing  templates  for  each  type  of  node  and  applying 
them  to  elements.  In  the  geometry  stage,  an  efficient  surface 
fitting  technique  is  developed  to  improve  surface  accuracy. 
The  constructed  T-splines  interpolate  every  boundary  node 
in  the  input  mesh,  with  C2 -continuity  everywhere  except  the 
local  region  around  irregular  nodes.  Finally,  Bezier  elements 
are  extracted  and  linear  independence  of  the  constructed  T- 
spline  is  studied  to  facilitate  isogeometric  analysis  [2,  10]. 

The  remainder  of  this  paper  is  organized  as  follows.  Sec¬ 
tion  2  reviews  T-splines  and  defines  rational  T-splines.  Sec¬ 
tion  3  explains  the  converting  algorithm  in  detail.  Section  4 
discusses  sharp  feature  preservation  and  surface  fitting.  Sec¬ 


tion  5  describes  a  Bezier  extraction  technique  to  facilitate 
isogeometric  analysis  and  studies  the  linear  independence 
of  T-splines.  Section  6  presents  results,  and  Section  7  draws 
conclusions. 

2  Rational  T-spline 

T-splines  [13]  are  generalized  from  NURBS  [8],  the  prevail¬ 
ing  industrial  standard  for  surface  modeling  in  Computer 
Aided  Design  (CAD)  and  Computer  Aided  Manufacturing 
(CAM).  Given  a  set  of  (m  +  1)  x  (n  +  1)  control  points  Qj 
(i  =  0,1,---  ,  m;  j  =  0, 1 ,  •  •  •  ,  n),  non-negative  weights  Wij  as¬ 
sociated  with  Qj,  degree  d  and  two  global  knot  vectors,  u  = 
[uq^  U\ ,  •  •  •  ,  Ufn+d ,  1]  and  V  [vo,  V\  ,  •  •  •  ,  Vn+d i  Vn+d+ 1]> 

the  NURBS  surface  is  defined  as 

m  n 

EEC;  WijNi4(u)Nj4(y) 

Sat(u,v)  =  ^^ - ,  (1) 

E  E  wijNi4(u)Nj4(v) 

i=0j=0 
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where  N^d{u)  and  Nj^ (v)  are  B-spline  basis  functions  de¬ 
fined  by  the  two  knot  vectors.  NURBS  provides  a  unified 
geometry  representation  of  standard  analytical  shapes  (e.g., 
conics)  and  free-form  shapes.  Additionally,  NURBS  are  in¬ 
variant  under  affine  as  well  as  perspective  transformations. 
However,  NURBS  have  two  drawbacks:  they  do  not  allow 
local  refinement  and  all  the  control  points  must  lie  topolog¬ 
ically  in  a  rectangular  grid.  To  overcome  these  two  draw¬ 
backs,  Sederberg  et  al  developed  T- splines  [13],  which  al¬ 
low  T-junctions  and  L-junctions  in  their  control  grid.  A  T- 
junction  terminates  a  row  or  column  of  control  points  in  the 
control  grid,  for  example  node  A  in  Figure  2.  An  L-junction 
terminates  a  row  and  a  column  of  control  points,  like  node 
B  in  Figure  2,  which  is  not  permitted  in  analysis-suitable 
T-splines  [7,  2,  11].  A  T-spline  surface  is  defined  by 


- ,  Men,  (2) 

I >A-M 

z=0 

where  w,-  is  the  weight  for  the  control  point  Q,  |)  = 
Nf(QN^(r\),Nf  and  /V-1 *  are  B-spline  basis  functions  defined 
by  two  local  knot  vectors,  ^  =  fei0,  t,n,  C/2,  C/3,  Cm]  and 
r|  i  =  [r|/o,  r|,|.r|,2,  r|/3,  r|,4]  when  degree  d  =  3,  and  £2  is  the 
local  domain1  of  the  T-spline  in  parameter  space. 


Fig.  2  One  local  region  of  a  T-mesh.  Node  A  is  a  T-junction  and  node 
B  is  an  L-junction. 


A  T-mesh  provides  the  connectivity  of  the  control  points 
and  a  knot  interval  is  assigned  to  each  edge  in  the  T-mesh  to 
indicate  the  parametric  length  of  that  edge.  The  local  knot 
vectors  for  each  node  are  inferred  from  the  T-mesh.  In  [12], 
Sederberg  et  al.  introduced  three  types  of  T-spline  spaces: 
standard,  semi-standard  and  non-standard,  based  on  whether 
the  basis  functions  or  the  weighted  basis  functions  can  pro¬ 
vide  a  partition  of  unity.  T-splines  provide  more  flexibili¬ 
ties  for  modeling.  However,  there  are  several  open  problems 

1  In  this  paper,  “domain”  refers  to  one  parametric  area  in  2D  or  one 

parametric  volume  in  3D  and  “patch”  refers  to  the  T-spline  surface  or 

solid  T-spline  defined  on  one  domain. 


which  limit  their  application,  such  as  how  to  characterize  T- 
mesh  configurations  for  a  standard,  semi- standard,  or  non¬ 
standard  T-spline,  and  how  to  calculate  the  weights  for  a 
semi-standard  T-spline. 

In  order  to  obtain  basis  functions  satisfying  a  partition 
of  unity,  we  choose  the  rational  basis  functions  to  construct 
T-splines.  The  rational  T-spline  surface  is  defined  as 


sM  =  ^ - -  Men, 

i= 0 


where 

*«M 


^(C)M) 

I”=0^(^(T1) 


(3) 


(4) 


is  the  newly  defined  rational  B-spline  basis  function,  N?  and 
N ^  are  B-spline  basis  functions  defined  by  the  local  knot 
vectors  at  node  Q,  ^  =  feo,  £/i,  ^2,  ^3,  and  rp  =  [rp0, 
rpi,  rp2,  Tp3,  rp4]  when  degree  d  =  3.  Similarly,  the  formula 
for  a  rational  solid  T-spline  is 


EwiCiftM  Q 

S&i i,C)  =  ^ - ,  MQen, 

1,0 

i= 0 


where 

*/MC) 


(5) 


(6) 


is  the  newly  defined  rational  B-spline  basis  function,  n}, 

r 

and  A-3  are  B-spline  basis  functions  defined  by  the  local  knot 
vectors  at  node  Q,  q,  =  [C/o,  C/i,  C/2,  Co,  Cm],  %  =  [%o,  T|n, 
B/2,  r|/3,  U/4]  and  Q  =  [£,-0,  C/i,  C/2,  C/3,  Cm]  when  degree 
J  =  3.  It  is  obvious  that  the  rational  B-spline  basis  functions 
automatically  satisfy  Y!j=o^i  =  1,  for  any  (^,r|)  in  2D  and 
(£,,rp  Q  in  3D.  In  this  way,  we  obtain  one  set  of  basis  func¬ 
tions  satisfying  a  partition  of  unity  even  for  non-standard 
T-splines,  and  successfully  avoid  the  difficulty  of  checking 
the  type  of  T-splines  from  the  T-mesh  configuration  and  cal¬ 
culating  the  weights  for  semi- standard  cases. 


3  Converting  Algorithm 

As  shown  in  Figure  3,  there  are  two  main  stages  in  con¬ 
verting  an  unstructured  quadrilateral  or  hexahedral  mesh  to 
a  T-spline  surface  or  solid  T-spline:  the  topology  stage  and 
the  geometry  stage.  We  take  the  input  mesh  as  the  initial 
T-mesh  and  the  topology  stage  aims  to  make  the  initial  T- 
mesh  gap-free  by  designing  templates  for  each  type  of  node 
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(a)  (b)  (c) 


Fig.  4  Regular  node,  partial  extraordinary  node  and  extraordinary  node  in  2D  and  3D.  Node  A  rendered  in  white  is  a  regular  node  (a,  d),  the 
magenta  node  is  a  partial  extraordinary  node  (e),  and  red  nodes  are  extraordinary  nodes  (b,  c,  f  and  g). 


Guad/HeK 

Mesh 

^ _ > 

Topology  stage 

f - ^ 

Gapfree 

T-mesh 

\ ^ 

Geometry  stage 

T-5pline 

c _ J 

*  Template 
Implementation 

[*  Surface  Fitting  J 

■  Sharp  Feature 

Preservation 

V, _ _ _ _ J 

Fig.  3  An  overview  of  the  algorithm  to  convert  an  unstructured  quadri¬ 
lateral  or  hexahedral  mesh  to  a  T-spline  surface  or  solid  T-spline. 

and  applying  them  to  elements.  Additional  nodes  and  edges 
are  inserted  to  preserve  sharp  features  in  the  input  mesh.  We 
assign  a  unit  knot  interval  to  each  edge  in  the  input  quadri¬ 
lateral  or  hexahedral  mesh.  All  the  new  edges  inserted  in 
the  topology  stage  have  either  zero  or  unit  parametric  edge 
length.  As  a  result,  a  quasi-uniform  T-mesh  is  constructed, 
which  contains  edges  with  only  zero  or  unit  knot  interval. 
Edges  with  zero  knot  interval  are  called  zero-length  edges. 
In  the  geometry  stage,  the  goal  is  to  minimize  the  error  be¬ 
tween  the  input  mesh  and  the  output  T-spline.  Then  a  T- 
spline  surface  or  solid  T-spline  is  constructed  based  on  the 
obtained  T-mesh.  In  addition,  Bezier  elements  are  extracted 
and  linear  independence  is  studied  in  order  to  facilitate  iso¬ 
geometric  analysis.  Here  are  some  definitions  which  will  be 
needed  in  the  following  algorithm  description. 

Definition  3.1.  A  pair  of  reflection  edges  are  two  adjacent 
edges  with  one  common  node  and  all  the  elements  sharing 
one  edge  are  topologically  symmetric  with  all  the  elements 
sharing  the  other  one,  with  respect  to  a  line  of  symmetry  in 
2D  or  a  plane  of  symmetry  in  3D.  The  line  of  symmetry  in 
2D  is  formed  by  all  the  adjacent  edges  of  the  shared  node 
except  for  these  two  edges.  The  plane  of  symmetry  in  3D 


is  formed  by  the  adjacent  quadrilaterals  of  the  shared  node 
which  do  not  contain  any  of  these  two  edges.  In  a  pair  of 
reflection  edges,  one  edge  is  also  called  the  reflection  edge 
of  the  other  one  about  the  shared  node  and  vice  versa. 

Let  us  take  node  A  and  its  adjacent  edges  AB  and  AC  in 
Figure  4(a)  as  an  example.  The  blue  line  DE  is  the  symme¬ 
try  line,  formed  by  the  adjacent  edges  of  A  except  for  AB 
and  AC.  The  two  quadrilaterals  adjacent  to  AB  and  the  two 
adjacent  to  AC  are  topologically  symmetric  with  respect  to 
the  symmetry  line  DE ,  hence  AC  is  the  reflection  edge  of  AB 
about  node  A  and  vice  versa.  For  node  A  and  its  two  adjacent 
edges  AB  and  AC  in  Figures  4(d-e),  the  blue  face  is  the  sym¬ 
metry  plane,  formed  by  the  quadrilateral  elements  adjacent 
to  node  A  but  not  containing  AB  or  AC.  The  hexahedral  ele¬ 
ments  adjacent  to  AB  and  all  the  elements  adjacent  to  AC  are 
topologically  symmetric  with  respect  to  the  blue  symmetry 
plane,  hence  AB  and  AC  are  a  pair  of  reflection  edges. 

Definition  3.2.  A  regular  node  is  a  node  about  which  each 
adjacent  edge  has  a  reflection  edge. 

In  Figure  4(a),  each  edge  adjacent  to  A  has  a  reflection 
edge.  By  definition,  node  A  is  a  regular  node  in  2D.  Node  A 
in  Figure  4(d)  is  a  regular  node  as  well,  because  each  edge 
adjacent  to  A  has  a  reflection  edge  about  A  in  3D.  A  regular 
node  always  has  a  valence  of  four  in  2D  and  a  valence  of 
eight  in  3D.  However,  a  node  with  a  valence  of  four  or  eight 
may  not  be  a  regular  node  in  general. 

Definition  3.3.  A  partial  extraordinary  node  is  an  irregular 
node  about  which  some  but  not  all  of  its  adjacent  edges  have 
reflection  edges. 

Let  us  take  node  A  in  Figure  4(e)  as  an  example,  its  two 
adjacent  edges,  AB  and  AC,  are  a  pair  of  reflection  edges. 
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However,  the  other  adjacent  edges  do  not  have  reflection 
edges.  Therefore,  by  definition  node  A  is  a  partial  extraor¬ 
dinary  node. 

Definition  3.4.  An  extraordinary  node  is  an  irregular  node 
about  which  none  of  its  adjacent  edges  has  a  reflection  edge. 

For  example,  node  A  in  Figure  4(b)  and  (c)  is  an  ex¬ 
traordinary  node.  In  2D,  an  extraordinary  node  has  a  valence 
other  than  four.  Figures  4(f)  and  (g)  show  two  extraordinary 
nodes  in  3D,  and  for  the  two  local  regions,  none  of  the  edges 
adjacent  to  A  has  a  reflection  edge.  In  2D,  we  can  use  the 
valence  number  to  identity  whether  a  node  is  regular  or  ex¬ 
traordinary.  However  in  3D,  situations  are  more  complicated 
and  there  is  no  such  direct  relationship.  In  3D,  a  regular  node 
has  a  valence  of  eight,  but  a  valence-eight  node  may  not  be 
regular.  For  example,  node  A  in  Figure  4(g)  has  a  valence  of 
eight  but  it  is  an  extraordinary  one. 

We  classify  all  the  nodes  in  the  input  mesh  into  three 
categories:  regular,  partial  extraordinary  and  extraordinary 
nodes,  and  treat  them  differently  during  template  design. 
Regular  nodes  do  not  introduce  gaps  in  the  T-spline  or  de¬ 
crease  the  surface  continuity.  If  all  the  nodes  in  the  input 
mesh  are  regular,  for  example  a  structured  mesh,  the  initial 
T-mesh  is  topologically  correct  and  gap-free.  Unlike  regular 
nodes,  partial  extraordinary  nodes  and  extraordinary  nodes 
need  to  be  handled  properly,  otherwise  they  may  introduce 
gaps  in  the  T-spline  model.  Hence,  a  rule  is  defined  for  the 
template  design:  for  each  partial  extraordinary  node  or  ex¬ 
traordinary  node,  the  template  should  ensure  the  constructed 
T-mesh  is  gap-free. 


(a)  (b)  (c) 

•  Original  extraordinary  node  —  Nonzero-length  edge 

3  Original  regular  or  extraordinary  node  — Zero-length  edge 

•  Inserted  node  Nonzero-area  parametric  domain- 

Zero-area  parametric  domain 

Fig.  5  The  general  templates  for  an  extraordinary  node  in  2D,  which 
are  used  to  design  templates  for  each  type  of  element  for  set  1  (a),  sets 
2  and  3  (b),  and  set  4  (c)  in  Table  1. 


3.1  Converting  a  Quadrilateral  Mesh  to  a  T-spline  Surface 

Figure  5  shows  three  general  templates  for  an  extraordi¬ 
nary  node  in  2D.  (a)  was  derived  using  T-NURCCs  (Non- 
Uniform  Rational  Catmull-Clark  Surfaces  with  T-junctions) 
[13],  (b)  and  (c)  are  two  simplified  templates  based  on  (a) 
with  fewer  newly  inserted  nodes  and  edges.  Based  on  these 
templates,  we  design  four  sets  of  templates  for  each  quadri¬ 
lateral  element  type,  see  Table  1 .  There  are  six  types  of  el¬ 


(a)  (b)  (c) 


Fig.  6  The  template  derivation  process  for  element  type  2  (set  1)  from 
the  template  in  Figure  5(a).  (a)  The  quadrilateral  element  with  two  ex¬ 
traordinary  nodes;  (b)  the  result  after  applying  the  template  to  the  ex¬ 
traordinary  node  A;  and  (c)  the  final  result  after  applying  the  template 
to  the  other  extraordinary  node  B. 


(a)  (b)  (c) 

Fig.  7  Three  possible  results  for  element  type  2  (sets  2  and  3)  after 
applying  the  template  in  Figure  5(b)  by  choosing  different  orientations. 


Fig.  8  The  template  derivation  process  for  element  type  2  (set  4)  from 
the  template  in  Figure  5(c).  (a)  The  result  after  applying  the  template 
to  the  extraordinary  nodes  A  and  B;  and  (b)  the  final  result  after  adding 
the  blue  edge  according  to  a  T-spline  rule. 

ements  in  the  initial  T-mesh  classified  by  the  number  of  ex¬ 
traordinary  nodes:  elements  with  none,  one,  two  (neighbor¬ 
ing  or  diagonal),  three  and  four  extraordinary  nodes. 

First,  let  us  take  set  1  as  an  example  to  see  how  to  de¬ 
rive  these  templates  in  Table  1  from  the  template  in  Figure 
5(a).  Figure  6  demonstrates  the  derivation  process  for  ele¬ 
ment  type  2.  We  first  apply  the  template  to  the  extraordinary 
node  A  to  obtain  the  result  in  Figure  6(b),  and  then  apply 
the  template  to  the  other  ordinary  node  B  to  get  the  final  re¬ 
sult.  The  other  templates  for  set  1  can  be  derived  in  the  same 
manner.  For  sets  2,  3  and  4,  we  can  proceed  similarly  to  de¬ 
rive  templates  for  each  type  of  element  using  the  general 
templates  in  Figure  5(b-c).  Note  that  the  template  in  Figure 
5(b)  has  two  possible  orientations.  We  always  choose  the 
one  which  introduces  a  minimum  number  of  newly  inserted 
nodes  and  edges.  Figure  7  shows  three  possible  results  for 
element  type  2  (sets  2  and  3)  after  applying  the  template  in 
Figure  5(b)  by  choosing  different  orientations.  Obviously, 
(a)  is  the  best  result  because  it  introduces  the  minimum  num¬ 
ber  of  new  nodes.  Figure  8  shows  the  template  derivation 
process  for  element  type  2  (set  4)  from  the  template  in  Fig- 
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Table  1.  Templates  for  six  quadrilateral  element  types. 


0  Original  regular  node  •  Inserted  node  —  Nonzero- length  edge 

#  Original  extraordinary  node  —  /ero-lenglh  edge 


Non  zero-area  parametric  domain 
Zero-area  parametric  domain 


ure  5(c).  Figure  8(a)  is  obtained  by  applying  the  template 
to  the  extraordinary  nodes  A  and  B.  The  blue  edge  in  Fig¬ 
ure  8(b)  is  added  according  to  a  T- spline  rule  [13]:  if  a  T- 
j  unction  or  L-j unction  on  one  edge  of  a  face  can  “legally” 
be  connected  to  a  T-junction  or  L-j  unction  on  an  opposing 
edge  of  this  face,  the  two  T-junctions  or  L-junctions  must  be 
connected  in  the  T-mesh.  “Legal”  here  means  that  the  sum 
of  knot  intervals  on  opposing  edges  of  any  face  must  always 
be  equal. 

Sets  1  and  2  were  given  in  [15],  and  set  2  was  simplified 
from  set  1  in  order  to  insert  fewer  nodes  and  get  better  con¬ 
tinuity.  In  [15],  the  designed  templates,  together  with  a  T- 
mesh  standardization  algorithm,  were  used  to  convert  an  un¬ 
structured  quadrilateral  mesh  to  a  standard  T- spline  surface. 
In  this  paper,  we  choose  the  rational  basis  functions  and  do 
not  need  to  get  a  standard  T-spline.  We  also  simplify  the  tem¬ 
plate  of  type  4  elements  in  set  2  to  obtain  set  3.  Sets  1,  2  and 
3  guarantee  a  gap-free  T-spline  surface  as  proved  in  [15]: 
for  a  T-mesh  without  L-junctions,  if  the  region  formed  by 
all  the  adjacent  elements  of  one  extraordinary  node  does  not 
contain  any  T-junctions  and  all  the  edges  in  it  have  zero  knot 
intervals,  the  local  region  around  this  extraordinary  node  is 
gap-free.  These  three  sets  work  for  any  T-mesh.  Since  we 


focus  on  quasi-uniform  T-meshes  in  this  paper,  we  simplify 
set  3  further  and  obtain  set  4,  which  is  specially  designed 
for  quasi-uniform  T-meshes  and  with  L-junctions  involved. 
Set  4  can  also  generate  gap-free  T-splines  as  stated  in  the 
following  lemma. 

Lemma  1  For  any  input  unstructured  quadrilateral  mesh, 
the  quasi-uniform  T-mesh  obtained  by  applying  template  set 
4  is  gap-free. 

Proof:  We  first  prove  this  lemma  for  a  valence-3  node  and 
then  extend  it  to  other  extraordinary  nodes.  As  shown  in  Ta¬ 
ble  2,  there  are  three  non-zero  domains  around  a  valence- 3 
extraordinary  node  O  and  we  first  prove  the  patch  defined  by 
£2°  is  continuous  or  gap-free  with  the  patch  defined  by  Q1. 
In  other  words,  we  need  to  prove  that  all  the  nodes  share  the 
same  basis  function  value  at  the  boundary  of  the  two  neigh¬ 
boring  domains  £2°  and  £2:.  In  Table  2,  the  top  row  pictures 
show  the  local  region  of  the  T-mesh  around  the  extraordi¬ 
nary  node  O.  The  second  and  third  rows  show  the  nodes 
around  the  extraordinary  node  with  non-zero  basis  function 
values  and  their  parametric  position  for  £2°  and  £^,  respec¬ 
tively.  Node  O  is  the  parametric  origin.  There  are  four  cases 
for  domains  £2°  and  £2*:  Nodes  L  and  J  are  both  extraordi- 
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Table  2.  Local  region  around  one  valence-3  extraordinary  node. 


Case  1 


Case  2 


Case  3 


Case  4 
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O  Regular  node  3  Regular  or  extraordinary  node 

6  Extraordinary  node 


Uml-lenglh  edge 
Zero- length  edge 


Nonzero-area  parametric  domain 
Zero-area  parametric  domain 


nary  nodes  (case  1);  either  L  or  /  is  an  extraordinary  node 
(case  2);  nodes  L,  J  and  K  are  all  regular  (case  3);  and  nodes 
L  and  J  are  regular  while  node  K  is  extraordinary  (case  4). 
The  type  of  nodes  G,  H ,  I  will  not  influence  the  continuity 
of  the  two  patches  defined  by  QP  and  £ll . 

In  case  1,  we  can  observe  that  except  for  node  O,  all 
the  other  nodes  with  non-zero  basis  function  values  in  £2° 
and  £2*  share  the  same  knot  vectors.  In  other  words,  they 
share  the  same  basis  function  in  these  two  domains.  In  addi¬ 
tion,  although  node  O  has  different  basis  functions  for  these 
two  domains,  it  shares  the  same  basis  function  value  at  the 
shared  boundary  of  the  two  domains  (^  =  0,T|  E  [—1,0]), 
because  N^  a 0(0)  =  (0)  =  1  and  A^o  =  a ,  when 

T|  E  [—1,0].  Therefore,  we  can  conclude  that  for  case  1,  the 
two  patches  defined  by  £2°  and  £ll  are  gap-free. 

In  case  2,  except  for  nodes  O  and  K ,  all  the  other  nodes 
with  non-zero  basis  function  values  in  £2°  and  share  the 
same  knot  vectors.  Similar  to  case  1,  node  O  shares  the  same 
basis  function  value  at  the  shared  boundary.  Node  K  has  zero 
basis  function  value  in  £2°,  and  although  it  has  non-zero  ba¬ 


sis  function  value  in  Q1,  the  function  value  is  zero  at  the 
shared  boundary.  Hence  again,  we  can  conclude  that  for  case 
2,  the  two  patches  are  gap-free. 

In  case  3,  nodes  O  and  K  do  not  share  the  same  knot 
vectors  in  £2°  and  £2X  and  the  other  nodes  share  the  same 
knot  vectors.  Again,  node  O  shares  the  same  basis  function 
value  at  the  shared  boundary.  Then  we  only  need  to  prove 
node  K  has  the  same  basis  function  or  the  same  basis  func¬ 
tion  value  at  the  shared  boundary.  Suppose  nodes  M  and  N 
are  the  other  two  adjacent  nodes  of  node  K ,  besides  L  and 
J.  There  are  three  possibilities  for  this  case  classified  by  the 
node  types  of  M  and  N :  both  of  them  are  regular;  both  of 
them  are  extraordinary;  one  is  regular  and  the  other  is  ex¬ 
traordinary. 

-  When  M  and  N  are  both  regular,  the  two  knot  vectors  of 
node  K  for  Q.0  are  qK.n(l  =  [-3,  -2,-1,  0, 1],  -r\K,a()  = 
[—  1 , 0, 1 , 2, 3]  and  the  two  knot  vectors  for  £ll  are  %k,cix  = 
[-1,0, 1,2, 3],ri^0l  =  [-1,0, 1,2, 3].  Since  vf  ^ (0)  = 
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Fig.  9  A  local  region  around  one  valence-5  extraordinary  node  O.  (b)  and  (c)  show  the  nodes  around  the  extraordinary  node  with  non-zero  basis 
function  values  and  their  parametric  position  for  domains  OP  and  . 


(0)  —  6’  we  ^ave  4^)N^(J\)  —  i(^) 

whenri  G  [-1,0]. 

-  When  M  and  N  are  both  extraordinary,  the  two  knot  vec¬ 
tors  of  node  K  for  £2°  are  ^k.q0  =  [— 2,  —  2,  — 1,0, 1], 
T|^^0  =  [— 1,0, 1,2,2]  and  the  two  knot  vectors  for  £ll 
are  =  [-1,0, 1,2,2],  r\ =  [-1,0, 1,2,2].  Since 

4a>(°)  =  4a(°)  =  5,  wehave^o^Xno^)  = 

4  a  (°)  4  a  (Tl)  when  n  e  [-1,0]. 


-  Suppose  M  is  extraordinary  and  N  is  regular.  The  two 
knot  vectors  of  node  K  for  £2°  are  ^k.q0  =  [—2,  — 2,  —  1, 
0,1],  T| k,q0  =  [—1,0, 1,2,3]  and  the  two  knot  vectors 


for^1  are^!  =  [-1,0, 1,2,3],  =  [-1,0, 1,2,2]. 

Since  ^  (0)  =  Qi  (0)  =  ^ ,  and  using  the  Oslo  knot 
insertion  algorithm  [3],  we  have  ^(rj)  =  ^  (r|), 

T|  G  [—1,0].  Hence  we  can  obtain  ^  q0(0)N^  Qq  (r|)  = 

4a  (°)  4  a  ^  whenrl  e  [-1,0]. 


In  conclusion,  for  all  the  three  possibilities,  node  K  always 
has  the  same  basis  function  values  at  the  shared  boundary 
of  the  two  adjacent  domains.  Therefore,  the  two  patches  de¬ 
fined  by  £1°  and  £ll  are  gap-free. 

In  case  4,  except  for  nodes  O ,  P ,  Q  and  R ,  all  the  other 
nodes  with  non-zero  basis  function  values  in  £2°  and  £ll 
share  the  same  knot  vectors.  In  addition,  node  O  has  the 
same  basis  function  value  at  the  shared  boundary  of  the  two 
domains  (^  =  0,r|  G  [—1,0]).  Node  P  has  non-zero  basis 
function  value  in  £2°,  but  the  function  value  is  zero  at  the 
shared  boundary  and  also  it  has  zero  basis  function  value  in 
f}1.  Node  R  has  the  same  situation  with  node  P.  The  knot 
vectors  of  node  Q  are  ^Q:a0  =  [-1, -1, -1,0, 1],  r| 0A)  = 
[-1,0, 1,1,1]  for  £2°,  and  £Ga  =  [—1,0, 1, 1, 1],  T|2a  = 
[-1,0, 1,1,1]  for  O1.  Since  4 A  (°)  =  4a  (°)  =  b  we 

have  44°)4a4)  =  ^  e  h1.0!- 

Therefore,  we  can  conclude  that  for  case  4,  the  two  patches 
defined  by  £2°  and  £ll  are  gap-free. 


Fig.  10  1  -ring  and  2-ring  neighborhoods  around  an  extraordinary  node 
O  in  the  T-mesh  after  applying  the  four  template  sets  in  Table  1 .  (a)  set 
1;  (b)  sets  2  and  3;  and  (c)  set  4. 


In  summary,  for  all  these  four  cases,  the  two  patches  de¬ 
fined  by  £1°  and  Q.1  share  the  same  curve  when  s  =  0  and 
t  G  [ — C\ ,  0] .  In  other  words,  the  two  patches  are  continu¬ 
ous  or  gap-free  across  the  shared  boundary.  Due  to  symme¬ 
try,  the  surface  is  gap-free  across  the  boundary  shared  by 
the  two  patches  defined  by  £2*  and  £22  and  likewise  for  Q? 
and  £2°.  Hence  the  surface  is  gap-free  for  the  local  region  of 
the  quasi-uniform  T-mesh  around  a  valence- 3  extraordinary 
node  obtained  by  applying  template  set  4.  Here,  we  utilize  a 
valence-3  node  in  the  proof,  but  the  proof  can  be  generalized 
to  other  valence  numbers.  Let  us  take  one  valence-5  node  O 
in  Figure  9(a)  as  an  example,  if  we  treat  nodes  K ,  L  and  J  in 
the  same  way  as  the  nodes  K ,  L  and  J  in  the  valence- 3  proof, 
the  proof  proceeds  in  the  same  way.  □ 

Discussion:  In  summary,  sets  1  and  2  were  designed  in 
[15]  to  obtain  a  standard  T-mesh,  sets  3  and  4  are  based  on 
the  rational  T-spline  basis  functions.  Sets  1,  2  and  3  work 
for  arbitrary  T-meshes  and  set  4  can  only  be  used  for  quasi¬ 
uniform  T-mesh.  Template  set  4  inserts  many  fewer  nodes 
and  zero-length  edges,  leading  to  a  better  surface  continuity. 
For  the  rational  T-spline,  we  can  use  all  four  sets  of  tem¬ 
plates,  and  we  do  not  need  the  standardization  step.  Thus 
we  successfully  avoid  the  propagation  in  [15],  and  the  re¬ 
sults  have  fewer  nodes  and  better  surface  continuity. 

Let  us  check  the  surface  continuity  around  an  extraordi¬ 
nary  node  after  applying  the  four  sets  of  templates.  We  de- 
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Fig.  11  A  comparison  of  using  template  sets  1-4.  (a)  The  rational  T-mesh  after  applying  set  1;  (b)  the  standard  T-mesh  after  applying  set  2  and 
standardization;  (c)  the  rational  T-mesh  after  applying  set  3;  and  (d)  the  rational  T-mesh  after  applying  set  4. 


fine  the  p-ring  neighborhood  around  an  extraordinary  node 
as  follows:  The  one-ring  neighborhood  consists  of  all  the  T- 
mesh  faces  adjacent  to  the  extraordinary  node.  The  two-ring 
neighborhood  consists  of  the  one-ring  neighborhood  plus  all 
the  T-mesh  faces  adjacent  to  faces  on  the  one-ring  neighbor¬ 
hood.  This  process  is  repeated  as  many  times  as  necessary 
to  construct  the  p-ring  neighborhood  [15].  Figure  10  shows 
1-ring  and  2-ring  neighborhoods  after  applying  the  four  tem¬ 
plate  sets.  Using  set  1,  the  surface  continuity  is  C°  at  the  knot 
coordinate  corresponding  to  the  extraordinary  point  O ,  and 
C°  across  the  shared  curve  of  the  nonzero  area  parametric 
domain  for  each  ring  until  the  4-ring  neighborhood.  Using 
set  2  or  3,  the  surface  continuity  is  again  C°  at  the  knot  co¬ 
ordinate  corresponding  to  O ,  C°  across  the  shared  curve  of 
the  nonzero  area  parametric  domain  for  each  ring  until  the 
3 -ring  neighborhood,  and  C1  outside  the  3 -ring  neighbor¬ 
hood  until  the  4-ring  neighborhood.  Using  set  4,  the  surface 
continuity  is  C°  at  the  knot  coordinate  corresponding  to  O, 
and  C°  across  the  shared  curve  of  the  nonzero  area  paramet¬ 
ric  domain  for  each  ring  until  the  3 -ring  neighborhood.  Be¬ 
yond  the  3 -ring  neighborhood,  the  surface  is  C2 -continuous. 
Therefore,  set  4  is  the  best  choice  for  us  due  to  its  better 
surface  continuity  and  fewer  nodes  introduced. 


Figure  11  shows  a  comparison  of  using  template  sets  1- 
4.  (a)  shows  the  rational  T-mesh  after  applying  set  1.  The 
blue  edges  are  added  according  to  a  T-spline  rule  as  dis¬ 
cussed  earlier.  After  applying  set  2  or  3,  it  is  also  possible 
that  some  additional  edges  must  be  inserted  according  to  this 
rule,  (b)  shows  a  standard  T-mesh  using  set  2,  and  the  green 
edges  are  inserted  during  the  standardization  step.  If  we  re¬ 
move  all  the  green  edges,  (b)  becomes  the  result  of  a  rational 
T-mesh  after  applying  set  2.  (c)  and  (d)  show  the  results  af¬ 
ter  applying  sets  3  and  4,  respectively.  It  is  obvious  that  (d) 
contains  many  fewer  nodes  and  zero-length  edges,  resulting 
in  better  surface  continuity. 


3.2  Converting  a  Hexahedral  Mesh  to  a  Solid  T-spline 

Hexahedral  meshes  contain  three  different  types  of  nodes: 
regular,  partial  extraordinary  and  extraordinary  nodes.  Here, 
we  first  design  general  templates  for  partial  extraordinary 
and  extraordinary  nodes,  and  then  apply  them  to  each  type  of 
element.  In  3D,  one  hexahedral  element  has  eight  nodes  and 
each  node  has  three  possible  types.  In  addition,  one  partial 
extraordinary  node  has  three  possible  orientations.  Hence, 
each  node  has  a  total  of  five  possibilities  and  if  we  classify 
the  elements  as  we  did  for  quadrilateral  meshes,  there  will 
be  5 8  =  390625  types  of  elements  without  considering  sym¬ 
metry  and  complementary.  Therefore,  instead  of  listing  all 
the  templates  for  each  type  of  element  as  in  Section  3.1  (Ta¬ 
ble  1),  here  we  use  several  examples  to  explain  how  to  apply 
the  general  templates  to  a  certain  element  type. 

Figure  12(a-c)  show  the  general  templates  for  a  partial 
extraordinary  node  in  hexahedral  meshes.  The  magenta  edge 
adjacent  to  the  partial  extraordinary  node  has  a  reflection 
edge  about  this  node.  Figure  12(a-c)  are  generalized  from 
Figure  5(a-c),  respectively.  Similar  to  Figure  5(b),  the  tem¬ 
plate  in  Figure  12(b)  also  has  two  possible  orientations.  Again, 
we  choose  the  one  which  introduces  fewer  new  nodes  and 
edges.  Obviously,  among  these  three  templates,  Figure  12(c) 
is  the  best  choice  due  to  its  minimum  number  of  newly  in¬ 
serted  nodes.  Basically,  we  apply  the  2D  templates  in  Fig¬ 
ure  5  to  the  face  perpendicular  to  the  edge  which  has  a  re¬ 
flection  edge  about  this  partial  extraordinary  node.  It  is  no 
doubt  that  each  node  has  one  unique  knot  vector  along  the 
£  direction.  We  have  proved  the  geometry  is  gap-free  on 
the  iso-parametric  ^,T|-plane  in  [15]  and  in  Lemma  1.  In 
other  words,  A^(^)A^(r|)  for  any  node  in  one  domain  has 
the  same  function  values  with  its  neighboring  domain  at 
the  shared  boundary.  Hence,  A^(^)A^(r|)A^(Q  also  has  the 
same  function  value  at  the  shared  boundary  of  two  neigh¬ 
boring  domains,  and  the  obtained  geometry  is  gap-free  after 
applying  templates  in  Figure  12(a-c). 

Figure  13  shows  an  example  of  one  element  with  two 
partial  extraordinary  nodes.  Nodes  A  and  B  are  partial  ex- 
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(a)  (b) 

•  Original  primal  extraordinary  node 

•  Original  extraordinary  node 

C»  Original  regular,  pan  ini  extraordinary 
or  extraordinary  node 


(C) 

•  Inserted  node 

I  Wonzttro-volume  parametric  domain 
'  Zero-volume  parametric  domain 


(d)  (e) 

—  Nonzero- length  edge 

—  Zero-length  edge 

—  Bdge  which  has  reflection  edge 


Fig.  12  The  general  templates  for  partial  extraordinary  nodes  (a-c)  and  extraordinary  nodes  (d-e)  of  hexahedral  meshes,  (a-c)  are  generalized  from 
templates  in  Figure  5(a-c),  respectively;  (d)  is  extended  from  the  template  in  Figure  5(a);  and  (e)  is  a  simplified  version  of  (d). 


(a) 


Fig.  13  One  element  with  two  partial  extraordinary  nodes,  A  and  B. 
(a)  The  input  hexahedral  element  in  which  edge  AB  has  one  reflection 
edge  about  A  and  edge  BC  has  one  reflection  edge  about  B;  and  (b)  the 
result  after  applying  the  template  in  Figure  12(c)  to  nodes  A  and  B. 


traordinary  nodes,  edge  AB  has  a  reflection  edge  about  A 
and  edge  BC  has  a  reflection  edge  about  B.  For  each  partial 
extraordinary  node,  we  apply  the  template  in  Figure  12(c) 
to  the  face  perpendicular  to  its  adjacent  edge  which  has  one 
reflection  edge  about  it.  The  blue  dash  edge  in  (b)  is  added 
according  to  the  T-spline  rule  as  discussed  earlier. 

Figure  12(d-e)  give  two  general  templates  for  an  extraor¬ 
dinary  node  in  hexahedral  meshes,  which  are  generalized 
from  Figure  5.  The  template  in  Figure  12(d)  has  one  single 
orientation  and  it  can  guarantee  the  obtained  T-mesh  gap- 
free  as  proved  in  the  following  lemma. 

Lemma  2  For  the  local  region  of  any  extraordinary  node 
in  the  input  unstructured  hexahedral  mesh,  the  T-mesh  ob¬ 
tained  by  applying  the  template  in  Figure  12(d)  is  gap-free. 


two  neighboring  domains,  say  £2°  and  £2!,  share  the  same 
face  at  the  shared  boundary.  In  other  words,  we  need  to 
check  if  we  have  5°(^,0,  Q  =  S^jOjQ,  where  \  E  [0,^] 
and  £  E  [—£4,0]. 

In  £2°  and  £2] ,  only  the  nodes  on  the  shared  face  OGDC 
have  non-zero  basis  function  values  at  the  shared  bound¬ 
ary  of  these  two  domains.  In  addition,  these  nodes  share 
the  same  knot  vector  along  the  r\  direction  except  for  O, 
which  is  T|^0  ^1  =  [— £3, 0,0,0,  £i].  The  knot  vectors  of  O  are 
floA)  =  [-23,0,0,0,0]  andr|oA  =  [0,0,0,0,£i].  We  notice 
that  N[-e3 , 0, 0, 0,  ef  (0)  =  N[-e3 , 0, 0, 0, 0]  (0)  =  N[  0, 0, 0, 
0,£i](0)  =  1.  Additionally,  all  the  nodes  with  non-zero  ba¬ 
sis  function  value  at  the  shared  boundary  of  £2°  and  £2X 
share  the  same  knot  vectors  along  the  £,  and  £  directions 
in  these  two  domains.  For  instance,  =  [0,0, 0,0,  ef\ 

and  a  A  =  [—24, 0, 0, 0, 0] .  In  other  words,  they  have  the 


same  basis  functions  for  these  two  parametric  directions, 
f  (5)A/£q1  (Q-  Thus,  we  can  conclude 

that  {%)  (QN^o (il)  =  A^i  ©  m?nI  (il), 


which  means  that  the  two  patches  defined  by  £2°  and  Q1 
share  the  same  boundary  face  when  C  E  [0,^2],  T|  =  0  and 
C  E  [—^4,0],  and  the  two  patches  are  continuous  or  gap-free 
across  the  shared  boundary.  Due  to  symmetry,  this  is  also 


true  for  the  pairs  £21-£22,  £2°-£22,  £21-£23,  £22-£23  and  £2°-£23. 


Hence  the  solid  T-spline  is  gap-free  for  the  local  domains 
around  the  extraordinary  node  O.  □ 


Proof:  To  prove  the  T-mesh  obtained  is  gap-free,  we  need  to 
find  out  whether  or  not  the  two  patches  defined  by  any  two 
neighboring  domains  around  one  extraordinary  node  share 
the  same  face.  Let  us  take  the  local  T-mesh  around  an  ex¬ 
traordinary  node  O  in  Figure  14  as  an  example,  where  node 
O  is  the  parametric  origin,  (b)  shows  the  T-mesh  after  ap¬ 
plying  the  template  in  Figure  12(d).  To  prove  the  solid  T- 
spline  patch  defined  by  the  local  domain  is  gap-free,  we 
need  to  check  whether  or  not  the  two  patches  defined  by 


Compared  to  Figure  12(d),  Figure  12(e)  is  a  simplified 
version  and  it  was  extended  from  Figure  5(b).  The  2D  tem¬ 
plate  in  Figure  5(b)  has  one  property:  for  one  domain  £2 
surrounding  one  extraordinary  node  O,  nodes  in  other  sur¬ 
rounding  domains  which  share  only  node  O  with  £2  always 
have  zero  basis  function  value  in  £2.  When  extending  it  to 
3D,  we  want  to  design  templates  that  inherit  this  property. 
Figure  12(e)  is  designed  based  on  this  principle.  Basically, 
for  one  extraordinary  node  in  a  certain  element  we  insert 
one  cube  whose  edges  all  have  zero  knot  intervals,  extend 
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Fig.  14  One  local  region  with  an  extraordinary  node  O.  (a)  shows  the  input  hexahedral  mesh  and  (b)  is  the  exploded  view  of  the  T-mesh  after 
applying  the  template  in  Figure  12(d).  Node  O  is  the  parametric  origin. 


the  cube  along  one  adjacent  edge  until  the  element  bound¬ 
ary  and  insert  one  plane  parallel  to  one  adjacent  plane.  For 
example,  for  the  extraordinary  node  A  in  Figure  12(e),  we 
insert  one  cube  AVWU-IJT S ,  extend  the  cube  along  one  ad¬ 
jacent  edge  AD  until  the  element  boundary,  and  insert  one 
plane  IKPR  parallel  to  its  adjacent  face  ABCD.  In  contrast 
with  Figure  12(d),  there  are  six  possible  variations  for  this 
template  by  choosing  different  adjacent  edges  for  extending 
the  cube  and  different  adjacent  faces  for  inserting  the  par¬ 
allel  plane.  We  always  choose  the  orientation  which  brings 
the  fewest  newly  inserted  nodes.  For  one  extraordinary  node 
in  a  certain  element,  if  one  of  its  three  adjacent  nodes  is  ex¬ 
traordinary  or  partial  extraordinary,  the  edge  containing  it 
will  be  chosen  to  extend  the  cube.  Once  we  fix  the  orien¬ 
tation  to  extend  the  cube,  there  will  be  two  options  left  for 
inserting  the  plane  along  the  two  faces  sharing  this  edge.  Af¬ 
ter  we  insert  the  cube,  we  only  need  to  insert  two  additional 
nodes  in  order  to  form  one  plane.  Here,  we  still  choose  the 
orientation  which  will  insert  fewer  nodes.  If  none  of  its  three 
adjacent  nodes  is  extraordinary  or  partial  extraordinary,  we 
check  the  face-diagonal  nodes.  If  one  of  its  face-diagonal 
nodes  is  extraordinary  or  partial  extraordinary,  we  choose 
the  face  containing  that  node  to  insert  the  plane. 

Figure  15  gives  one  example  demonstrating  how  to  ap¬ 
ply  the  two  templates  in  Figure  12(d-e)  to  one  element  with 
two  extraordinary  nodes,  A  and  B.  (b)  shows  the  result  af¬ 
ter  applying  the  template  in  Figure  12(d)  for  node  B ,  and 
then  using  the  template  again  for  node  A  to  obtain  the  re¬ 
sult  shown  in  (c).  Since  we  have  one  unique  orientation,  this 
process  is  quite  straightforward,  (d)  shows  the  result  after 
applying  the  template  in  Figure  12(e)  for  node  B.  It  is  ob¬ 


tained  by  first  inserting  one  cube  around  B ,  then  extending 
the  cube  along  its  adjacent  edge  AB.  Since  node  A  is  ex¬ 
traordinary,  we  choose  edge  AB  to  extend  here.  After  that, 
we  have  two  options  to  choose  from  to  get  one  adjacent  face 
for  inserting  one  parallel  plane,  the  front  face  or  the  bottom 
face.  In  this  case,  both  of  these  options  will  end  up  with  in¬ 
serting  two  additional  nodes,  hence  we  can  choose  either  one 
of  them.  Here  the  bottom  face  is  chosen  and  plane  CDEF  is 
inserted.  Then  the  result  shown  in  (d)  is  obtained.  Similarly, 
since  node  A  already  has  one  plane  and  one  extended  cube, 
we  only  need  to  insert  one  cube  around  it.  Then  the  result  in 
(e)  is  obtained. 

Figure  16  gives  another  example  demonstrating  how  to 
apply  the  templates  in  Figure  12(c-e)  for  one  element  with 
one  extraordinary  node  B  and  one  partial  extraordinary  node 
C.  The  magenta  edge  has  one  reflection  edge  about  node  C. 
(b)  is  the  result  after  applying  the  template  in  Figure  12(d) 
for  extraordinary  node  B.  Then  we  apply  the  template  in  Fig¬ 
ure  12(c)  for  node  C  to  get  (c).  (d)  is  the  result  after  applying 
the  template  in  Figure  12(e)  for  node  B.  Here  we  choose  the 
front  face  to  insert  one  plane  because  node  C  lies  on  this 
front  face  diagonal  with  node  B.  (e)  shows  the  final  result 
after  applying  the  template  in  Figure  12(c)  for  node  C. 

Discussion:  In  summary,  we  design  three  general  tem¬ 
plates  for  a  partial  extraordinary  node  and  two  general  tem¬ 
plates  for  an  extraordinary  node.  Among  the  three  templates 
for  a  partial  extraordinary  node,  Figure  12(c)  introduces  many 
fewer  newly  inserted  nodes.  For  extraordinary  nodes,  we 
prove  that  using  the  general  template  in  Figure  12(d)  will 
produce  a  gap-free  solid  T- spline.  Figure  12(e)  is  simpli¬ 
fied  from  Figure  12(d),  and  due  to  its  complexity  we  do  not 
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(a)  (b)  (c)  (d)  (e) 


Fig.  15  An  example  demonstrating  how  to  apply  the  templates  in  Figure  12(d-e)  to  one  element  with  two  extraordinary  nodes,  A  and  B.  (a)  The 
hexahedral  element  with  two  extraordinary  nodes;  (b)  the  result  after  applying  the  template  in  Figure  12(d)  for  node  B\  and  (c)  the  final  result  after 
applying  the  template  in  Figure  12(d)  for  node  A.  Similarly,  (d-e)  shows  the  two  results  after  applying  the  template  in  Figure  12(e)  first  for  node  B 
and  then  for  node  A. 


Fig.  16  An  example  demonstrating  how  to  apply  the  templates  in  Figure  12(c-e)  to  one  element  with  one  extraordinary  node  B  and  one  partial 
extraordinary  node  C.  (a)  The  hexahedral  element  in  which  the  magenta  edge  has  one  reflection  edge  about  node  C;  (b)  the  result  after  applying 
the  template  in  Figure  12(d)  for  node  B;  and  (c)  the  final  result  after  applying  the  template  in  Figure  12(c)  for  node  C.  Similarly,  (d-e)  show  the 
two  results  after  applying  the  template  in  Figure  12(e)  first  for  node  B  and  then  the  template  in  Figure  12(c)  for  node  C. 


prove  that  it  can  always  guarantee  a  gap-free  solid  T-spline, 
although  our  experience  indicates  it  does. 

In  3D,  things  are  much  more  complicated  and  here  we 
only  discuss  the  continuity  of  the  neighboring  domains  ad¬ 
jacent  to  one  partial  extraordinary  or  extraordinary  node. 
In  Figure  12(a-c),  the  patch  defined  by  each  domain  is  C2- 
continuous  with  the  patch  defined  by  the  non-zero  domain 
sharing  the  bottom  face  with  it,  and  C° -continuous  with  the 
patches  defined  by  the  domains  sharing  the  front  or  the  right 
face.  In  Figure  12(d-e),  the  patch  defined  by  each  domain  is 
C° -continuous  with  the  patches  defined  by  all  the  non-zero 
domains  sharing  one  face  with  it. 

4  Sharp  Feature  Preservation  and  Surface  Fitting 

4.1  Sharp  feature  preservation 

For  T-spline  surfaces,  how  to  preserve  sharp  feature  was  de¬ 
scribed  in  [15].  The  main  idea  is  to  use  repeating  knots  to  de¬ 
crease  the  local  boundary  surface  continuity  to  C°,  and  zero- 
length  edges  are  inserted  across  sharp  edges  and  around  sharp 
corners  to  preserve  these  sharp  features.  For  each  sharp  edge, 
we  duplicate  it  for  each  face  and  all  the  transverse  edges 
have  zero  knot  interval.  Each  sharp  corner  is  treated  as  one 
extraordinary  node. 

For  solid  T- splines,  the  input  hexahedral  meshes  in  this 
paper  were  generated  using  an  octree-based  isocontouring 


Fig.  17  Sharp  feature  preservation  for  solid  T-spline.  (a)  Sharp  edge 
preservation  (the  blue  edge)  and  (b)  sharp  corner  preservation  (the  ma¬ 
genta  point). 

algorithm  together  with  a  pillowing  technique  [9].  Pillowing 
is  a  sheet  insertion  method  that  refines  the  mesh  boundary. 
After  pillowing  each  element  has  at  most  one  face  lying  on 
the  boundary.  For  each  boundary  element  we  first  insert  one 
face  parallel  to  the  boundary  face  and  the  edges  connecting 
them  have  zero  knot  intervals.  In  this  way,  only  the  bound¬ 
ary  nodes  have  non-zero  basis  function  value  on  the  solid 
T-spline  boundary.  For  example  in  Figure  17,  the  pink  faces 
are  boundary  faces  and  the  two  faces  GHKL  and  HUK  are 
newly  inserted  faces.  The  edges  shown  in  red  between  the 
boundary  faces  and  the  inserted  faces  have  zero  knot  inter¬ 
vals. 

To  preserve  sharp  edges  in  solid  T- splines,  we  insert  one 
edge  parallel  to  the  sharp  edge  for  each  adjacent  bound¬ 
ary  face.  In  contrast  with  sharp  feature  preservation  for  2D 
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Fig.  18  Surface  fitting  for  a  solid  T-spline  model,  (a)  The  input  hexahedral  mesh;  (b)  the  constructed  solid  T-spline  and  T-mesh  before  surface 
fitting;  and  (c)  the  constructed  solid  T-spline  and  T-mesh  after  surface  fitting. 


T-spline  surfaces,  we  need  to  repeat  inserting  zero-length 
edges  for  the  newly  inserted  faces.  For  example,  in  Figure 
17(a),  the  blue  edge  is  one  sharp  edge.  We  first  insert  two 
edges  on  the  boundary  faces  ABEF  and  BCDE.  This  step 
guarantees  there  are  two  zero-length  edges  across  the  sharp 
edge.  In  other  words,  there  are  three  duplicated  knots  and 
the  surface  continuity  is  C°  across  the  sharp  edge.  After  that, 
two  additional  edges  are  inserted  on  the  newly  inserted  faces 
GHKL  and  HIJK.  The  purpose  of  this  step  is  to  make  the 
boundary  face  and  the  layer  right  below  it  have  the  same 
topology. 

To  preserve  sharp  corners  in  solid  T- splines,  the  sharp 
corner  node  is  treated  as  one  partial  extraordinary  node.  Its 
adjacent  edge  containing  one  interior  node  is  treated  in  the 
same  way  as  when  it  has  one  reflection  edge.  The  templates 
in  Figure  12(a-c)  can  be  used  to  handle  the  sharp  comers. 
Again,  the  operation  done  for  the  boundary  face  needs  to  be 
repeated  for  the  newly  inserted  faces.  For  example,  in  Figure 
17(b)  the  magenta  node  represents  one  sharp  corner  and  the 
template  in  Figure  12(c)  is  applied.  Then  we  apply  the  same 
inserting  operation  to  the  pillowed  faces  GHKL  and  HIJK. 
In  this  way,  the  sharp  comer  is  preserved. 

4.2  Surface  fitting 

In  the  surface  fitting  step,  we  aim  to  relocate  T-mesh  nodes 
so  that  the  output  T-spline  interpolates  all  the  boundary  nodes 
in  the  input  mesh.  The  surface  fitting  algorithm  given  here 
for  2D  is  the  same  as  the  method  given  in  [15].  We  take  ad¬ 
vantage  of  the  local  property  of  T-splines  and  fit  each  Bezier 
element  to  interpolate  the  four  nodes  of  the  corresponding 
quadrilateral  element.  For  each  nonzero  domain  of  the  T- 
spline,  we  set  the  positions  of  the  four  nodes,  which  come 
from  the  input  mesh  boundary  and  whose  parametric  coor¬ 
dinates  correspond  to  the  four  corners  of  the  domain,  as  un¬ 
known  variables  and  fix  all  the  other  nodes.  Then  we  use  the 
interpolation  condition  to  calculate  the  new  position  of  the 


unknown  variables.  After  that,  all  the  nodes  in  the  T-mesh 
which  have  the  same  parametric  position  with  the  unknown 
variables  are  updated.  We  loop  over  each  non-zero  domain 
and  iterate  until  the  interpolation  error  falls  below  a  given 
tolerance. 

For  solid  T-splines,  the  boundary  nodes  are  relocated  in 
the  same  way  as  for  a  T-spline  surface.  After  each  reposi¬ 
tion  of  the  boundary  nodes,  we  also  need  to  relocate  the 
neighboring  interior  nodes  to  avoid  tangling  the  local  con¬ 
trol  mesh  and  to  improve  the  quality  of  the  Bezier  elements. 
Here,  we  borrow  some  ideas  from  the  subject  of  mesh  qual¬ 
ity  improvement.  The  interior  nodes  which  come  from  the 
input  mesh  are  relocated  using  smoothing  and  optimization, 
and  the  newly  inserted  interior  nodes  are  relocated  corre¬ 
spondingly.  In  smoothing,  each  interior  node  is  moved  to¬ 
wards  the  mass  center  using  its  neighboring  elements.  If 
the  smoothing  operation  cannot  improve  the  quality,  we  use 
optimization  for  this  node.  The  optimization  method  loops 
over  all  the  adjacent  hexahedral  elements  to  compute  their 
Jacobians,  and  then  the  element  with  the  worst  Jacobian  is 
found  and  improved  using  a  conjugate  gradient  method,  in 
which  the  objective  function  is  the  Jacobian  at  that  vertex 
[16].  The  Jacobian  is  used  as  a  metric  to  measure  the  mesh 
quality.  Note  that  relocating  control  points  to  improve  the 
quality  of  Bezier  elements  is  not  straightforward.  Here  we 
use  a  mesh  metric  to  relocate  the  control  points,  since  intu¬ 
itively  a  control  mesh  with  better  mesh  quality  will  yield  a 
T-spline  with  better  Bezier  elements.  Figure  18  shows  one 
surface  fitting  example  for  a  solid  T-spline,  in  which  the 
boundary  surface  of  the  solid  T-spline  after  surface  fitting 
in  Figure  18(c)  interpolates  all  the  boundary  nodes  in  the 
input  mesh. 

5  Bezier  Extraction  and  Linear  Independence 

To  facilitate  isogeometric  analysis,  Bezier  elements  are  ex¬ 
tracted  from  the  constructed  T-spline  surface  or  solid  T-spline 
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[2,  10].  For  each  nonzero  parametric  domain,  we  determine 
the  nodes  with  nonzero  basis  function  values  in  this  domain 
and  then  calculate  the  transformation  matrix  Me  between  the 
T- spline  basis  functions  and  the  Bezier  basis  functions.  For 
a  T-spline  surface,  we  have 


the  finite  element  method  to  get  a  global  transformation  ma¬ 
trix.  In  other  words,  we  calculate  the  global  transformation 
matrix,  K ,  for  the  whole  model  from  T-spline  to  Bezier  ele¬ 
ments.  Then  we  can  obtain 


Bet  =  MeB 


e 

bi 


(7)  B,  =  KBb, 


(13) 


where 


(8) 


is  the  vector  formed  by  the  T-spline  basis  functions  with 
nonzero  function  values, 


"MO,  0,0,0, 1]  (^^[0,0, 0,0,  l](r|) " 

M0, 0,0, 1,1]  (^[0,0, 0,0,1]  (11) 

Mo,  o,  i ,  i ,  i]  (S)Mo,  o,  o,  o,  i]  (r|) 


Mo,o,i,i,i](5)Mo,i,m,i](ti) 


(9) 


is  the  vector  formed  by  the  Bezier  basis  functions,  and  rf  is 
the  number  of  nodes  with  nonzero  basis  function  values  in 
this  domain.  Me  can  be  calculated  using  the  Oslo  knot  inser¬ 
tion  algorithm  [3].  Similarly,  for  a  solid  T-spline  we  have 

Bet=MeBeb,  (10) 


where 


where 


Bt  = 


N 11  N ^  N 7V^  1 

2VV0  ’7VV1  ^  >yV/i-2 


(14) 


is  the  vector  formed  by  all  the  T-spline  basis  functions2  of 
the  T-mesh,  and 


Bb 


[B°[0],Bg[l],...,B»[15],Bi[0],Bi[l],...,Bj[15], 

(15) 


is  formed  by  all  the  Bezier  basis  functions  of  the  model,  n  is 
the  number  of  nodes  in  the  T-mesh  and  m  is  the  number  of 
non-zero  domains  or  Bezier  elements  for  the  whole  model. 
The  matrix  K  is  the  global  transformation  matrix  from  T- 
spline  to  Bezier,  with  a  size  of  n  x  16m  in  2D.  In  3D, 


Bt  = 


1  T 


(16) 


and 


Bet  = 


and 


isfi N ^ N **  •  •  •  N ^  N 11  N ^  ]T 

,  mV-FV-FV-1  , 


(ID 


B6«[Bg[0],Bg[l],...,Bg[63],Bj[0],Bi[l],...,Bi[63], 


(17) 


'M0,  o,  0, 0,  m)N[0, 0, 0, 0,  l](r|)7V[0, 0, 0,0,1]  (Q " 
M0, 0,0, 1, 1]  ©M0, 0,0,0,  l](r|)7V[0, 0,0,0, 1](Q 
M0,0, 1, 1, 1]©M0, 0,0,0,  1](ti)A1[0, 0,0,0, 1](Q 


M0.0,  i,  i,  1](S)M0, 1, 1, 1, 1](T1)M0, 1, 1, 1,  i](Q 
_M0,  i,  1, 1, 1](S)M0, 1, 1, 1,  i](ri)M0, 1, 1, 1,  i](Q. 

(12) 


The  Bezier  extraction  technique  can  also  be  used  to  study 
the  linear  independence  of  a  given  T-mesh.  In  [7],  the  linear 
independence  of  a  T-spline  model  is  determined  by  comput¬ 
ing  the  nullity  of  the  T- spline- to-NURBS  transform  matrix. 
However,  this  method  is  not  suitable  for  a  T-spline  with  ex¬ 
traordinary  nodes,  since  converting  a  T-spline  with  extraor¬ 
dinary  nodes  to  NURBS  will  end  up  with  multiple  NURBS 
patches.  The  T-meshes  in  this  paper  have  a  large  percent¬ 
age  of  extraordinary  nodes.  In  addition,  we  need  to  check 
the  linear  independence  for  3D  solid  T- splines.  Here,  in  or¬ 
der  to  study  the  linear  independence  of  T-splines  we  assem¬ 
ble  all  transformation  matrices  of  the  non-zero  domains  or 
Bezier  elements  in  a  similar  way  as  for  matrix  assembly  in 


The  matrix  K  is  the  global  transformation  matrix  from  solid 
T-spline  to  solid  Bezier  and  the  size  of  K  is  n  x  64m.  Given 
one  T-mesh,  all  the  T-spline  basis  functions  form  a  linear 
space.  If  all  these  functions  are  linearly  independent,  they 
form  a  basis  of  the  space  with  dimension  n.  The  following 
Lemma  3  provides  a  necessary  and  sufficient  condition  for 
linear  independence  of  a  T-spline. 

Lemma  3  (Necessary  and  Sufficient  Condition  for  Lin¬ 
ear  Independence)  A  necessary  and  sufficient  condition  for 
a  T-spline  model  to  be  linearly  independent  is  that  the  global 
transformation  matrix  from  T-spline  to  Bezier  has  full  rank. 

Proof:  By  definition,  the  T-spline  basis  functions  are  linearly 
independent  if  and  only  if  there  do  not  exist  scalars,  a  = 
[oco,  •  *  *  ,  oc n]T ,  not  all  zero,  such  that 

aT  B,  =  0.  (18) 

2  The  term  “basis  function”  implies  linear  independence  and  it  is 
more  appropriate  to  use  the  terminology  “blending  function”.  How¬ 
ever,  since  the  majority  of  T-meshes  are  linearly  independent,  the  term 
“basis  function”  is  used  here  for  simplicity. 
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Table  3.  Statistics  of  all  the  tested  models 


Model 

Input  mesh 
(vertices,  elements) 

(PENs,  ENs) 

Number  of 
T-mesh  nodes 

Number  of 
Bezier  elements 

Time 

(s) 

Genus-3  model 

(3,068,  3,072) 

(0,  4) 

3,132 

3,072 

0.24 

Head 

(2,909,  2,908) 

(0,  1,222) 

12,677 

2,908 

3.5 

Quad 

Hand 

(6,070,  6,068) 

(0,  2,527) 

26,270 

6,068 

8.1 

Ribosome  30S 

(13,217,  13,215) 

(0,  7,217) 

70,937 

13,215 

25.2 

Ribosome  50S 

(16,537,  16,537) 

(0,  9,039) 

88,849 

16,537 

27.7 

Statue 

(14,095,  12,313) 

(2,738,  1,663) 

114,672 

12,313 

57.9 

Hex 

Bearing 

(12,184,  10,215) 

(2,535,  756) 

71,339 

10,215 

39.8 

Gear 

(23,642,  19,438) 

(6,149,  2,021) 

171,827 

19,438 

83.1 

Eddy  plate 

(28,747,  24,269) 

(8,309,  2,906) 

239,660 

24,269 

111.2 

Note:  PEN  indicates  partial  extraordinary  node  and  EN  indicates  extraordinary  node. 


In  other  words, 

aTKBb  =  0.  (19) 

Since  the  Bezier  basis  functions  are  linearly  independent, 
the  necessary  and  sufficient  condition  for  linear  dependence 
of  the  T-spline  basis  functions  becomes 

a  tK  =  Kt  a  =  0,  (20) 

for  a  /  0.  This  will  only  happen  when  K  does  not  have  full 
rank.  □ 

Discussion:  Lemma  3  utilizes  the  transformation  matrix 
from  T-spline  to  Bezier  elements  to  study  the  linear  inde¬ 
pendence  of  a  given  T-spline.  This  technique  is  more  gen¬ 
eral  compared  with  the  linear  independence  study  using  the 
T-spline-to-NURBS  transformation  matrix.  It  works  for  T- 
meshes  with  extraordinary  nodes.  However,  generally  the 
size  of  matrix  K  is  very  large,  and  checking  its  rank  can  be 
very  time-consuming. 

6  Results 

We  used  template  set  4  in  Table  1  and  applied  the  converting 
algorithm  to  several  unstructured  quadrilateral  meshes  (Fig¬ 
ures  19-23).  If  the  input  mesh  contains  few  extraordinary 
nodes,  for  example  Figure  19,  the  output  T-spline  surface 
will  be  very  smooth.  We  also  applied  our  converting  algo¬ 
rithm  to  hexahedral  meshes  of  one  statue  model  and  several 
CAD  assembly  models  with  sharp  features  (Figures  24,  1, 
25  and  26).  For  these  3D  hexahedral  meshes,  the  templates 
in  Figure  12(c)  and  (d)  are  utilized  for  partial  extraordinary 
nodes  and  extraordinary  nodes,  respectively.  It  is  apparent 
that  the  constructed  solid  T-splines  preserve  all  the  sharp 
features  in  the  input  model.  In  Figures  1,  25  and  26,  differ¬ 
ent  colors  represent  different  components  of  the  assembly 
model.  These  components  have  conformal  boundaries  in  the 
output  solid  T-spline.  The  converting  algorithm  is  efficient 
and  all  the  results  were  computed  on  a  PC  equipped  with 
an  Intel  Q6600  (4  cores,  2.4GHz)  processor  and  4GB  main 
memory  (DDR2,  800MHz). 


Statistics  for  all  the  tested  models  are  shown  in  Table  3. 
From  the  table,  we  can  notice  that  the  number  of  Bezier  el¬ 
ements  is  the  same  as  the  number  of  elements  in  the  input 
mesh.  The  reason  for  this  is  all  the  edges  we  inserted  have 
zero  knot  interval.  There  is  a  subtlety  here  that  is  difficult 
to  see.  The  input  mesh  consists  of  bilinear  quadrilateral  ele¬ 
ments  or  trilinear  hexahedral  elements.  The  output  T-spline 
surface  is  bicubic  and  C2 -continuous  except  in  the  vicinity 
of  extraordinary  nodes  and  the  solid  T-spline  is  tricubic  and 
C2 -continuous  except  in  the  vicinity  of  partial  extraordinary 
and  extraordinary  nodes.  The  Bezier  elements  are  embed¬ 
ded  in  the  T-spline.  For  all  the  tested  models,  the  output  T- 
spline  surface  or  solid  boundary  interpolates  all  the  nodes  in 
the  input  quadrilateral  mesh  or  hexahedral  mesh  boundary. 
The  input  hexahedral  meshes  of  the  three  CAD  models  were 
generated  from  NURBS  boundary  representations  in  [9].  Al¬ 
though  the  output  solid  T-splines  interpolate  all  the  nodes  on 
the  mesh  boundary,  there  are  differences  compared  with  the 
original  NURBS  model. 


7  Conclusions 

We  have  developed  a  novel  algorithm  for  converting  any  un¬ 
structured  quadrilateral  or  hexahedral  mesh  to  a  T-spline  sur¬ 
face  or  solid  T-spline,  respectively.  The  T-spline  definition  is 
based  on  rational  T-spline  basis  functions,  with  the  partition 
of  unity  property.  The  converting  algorithm  has  two  main 
stages:  the  topology  stage  and  the  geometry  stage.  In  the 
topology  stage,  we  design  templates  for  each  type  of  node 
and  element  in  order  to  get  a  gap-free  T-spline.  Sharp  fea¬ 
tures  are  preserved  automatically  in  this  stage.  In  the  geom¬ 
etry  stage,  an  efficient  surface  fitting  technique  is  developed 
to  improve  the  surface  accuracy.  Finally,  a  Bezier  extraction 
technique  is  introduced  and  linear  independence  of  the  con¬ 
structed  T-spline  is  studied  to  facilitate  T-spline  based  iso¬ 
geometric  analysis.  As  part  of  future  work,  we  plan  to  study 
linear  independence  directly  from  the  T-mesh  configuration 
and  construct  solid  T-splines  directly  from  the  NURBS  bound¬ 
ary  representation. 
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Fig.  19  One  genus-3  model,  (a)  The  input  unstructured  quadrilateral  mesh;  (b)  the  constructed  T-spline  surface  and  T-mesh;  (c)  the  extracted 
Bezier  elements;  and  (d)  the  T-spline  surface,  (e)  to  (h)  show  details. 


Fig.  20  Human  head,  (a)  The  input  unstructured  quadrilateral  mesh;  (b)  the  constructed  T-spline  surface  and  T-mesh;  and  (c)  the  extracted  Bezier 
elements. 
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Fig.  21  Human  hand,  (a)  The  input  unstructured  quadrilateral  mesh;  (b)  the  constructed  T-spline  surface  and  T-mesh;  and  (c)  the  extracted  Bezier 
elements. 


Fig.  22  Ribosome  30S.  (a)  The  input  unstructured  quadrilateral  mesh;  (b)  the  constructed  T-spline  surface  and  T-mesh;  and  (c)  the  extracted  Bezier 
elements. 
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(a)  (b)  (c) 


Fig.  23  Ribosome  50S.  (a)  The  input  unstructured  quadrilateral  mesh;  (b)  the  constructed  T-spline  surface  and  T-mesh;  and  (c)  the  extracted  Bezier 
elements. 


Fig.  24  One  statue  model,  (a)  The  input  unstructured  hexahedral  mesh;  (b)  the  constructed  solid  T-spline  and  T-mesh;  and  (c)  the  extracted  solid 
Bezier  elements  with  some  elements  removed  to  show  the  interior  mesh. 
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